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Abstract 

The mathematical structure of the six dimensional physical phase spaces of the non-diagonal 
Bianchi IX model is analyzed in the neighborhood of the cosmological singularity. Critical points 
of the Hamiltonian equations appearing at infinities are of the nonhyperbolic type. Specific trans- 
formations of the phase space, including projection into finite region, do not change this type of 
criticality which is difficult for investigation by standard analytical tools. The nonhyperbolicity 
seems to be a generic feature of considered singular dynamics. The information that can be ob- 
tained from the linearized vector field is inconclusive. Making use of the physical Dirac observables 
as the phase space coordinates lowers substantially the dimensionality of the dynamics arena. Here, 
Qi using commonly known methods for studying the dynamics may turn out to be quite satisfactory. 
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I. INTRODUCTION 



The available data of observational cosmology indicate that the Universe emerged from a 
state with extremely high energy densities of matter fields. Theoretical cosmology, particu- 
larly the Belinskii-Khalatnikov-Lifshitz (BKL) scenario (general solution of general relativ- 
ity), predicts the existence of the initial cosmological singularity with diverging gravitational 
and matter fields invariants [H, 01 • It has been proved that such symmetries of spacetime like 
isotropy and homogeneity are dynamically unstable in the evolution towards the singularity 
[l|, 0]. An existence of the dynamically stable general solution to the Einstein equations with 
a cosmological singularity means that this classical theory is incomplete^. It is expected 
that finding the singularity free quantum cosmology would help in the construction of a the- 
ory unifying gravitation and quantum physics. The standard model of cosmology is unable 
to explain, at the fundamental level, the nature of dark energy and the origin of cosmic 
inflation. It is expected that quantum BKL theory may bring solution to these problems. 

The best prototype for the BKL scenario is the non-diagonal Bianchi IX model j2l, [sf. 
We expect that obtaining quantum Bianchi IX model will enable quantization of the BKL 
theory with the general cosmological singularity. Finding the non-singular quantum BKL 
theory would mean solving the generic cosmological singularity problem. Such a quantum 
theory could be used as a realistic model of the very early Universe. Obtaining quantum 
BKL theory may have fundamental importance for physics and cosmology. 

In this paper we examine the mathematical structure of the physical phase spaces of 
the non-diagonal Bianchi IX model near the cosmological singularity. In our analyses we 
apply the dynamical systems methods. We try to find some convenient variables in the 
phase space such that the critical points of the vector field specifying the dynamics are 
possibly of hyperbolic type. Quantum versions of these variables will be used as building 
blocks for finding quantum compound observables. The latter are supposed to be measurable 
observables. 

In Sec. II we identify the Hamiltonian of our system, which turns out to be the dynam- 
ical constraint, and we show that Hamilton's equations are equivalent to the Lagrangian 
equations. In Sec III we analyze the structure of the phase space and identify the subspace 
of critical points. It turns out to be of nonhyperbolic type. In the rest of the paper we con- 
sider two mappings of the phase space to find the one leading to hyperbolic type of critical 
points. The McGehee type transformation considered in Sec IV leads to another system with 
nonhyperbolic critical subspace. The commonly used transformation of the Poincare type, 
presented in Sec V, leads to the nonhyperbolic critical subspace as well. Finally, in Sec. VI 
we propose to use the Dirac observables to paramertize the phase space. We conclude in the 
last section. The Appendix presents the dynamics of the Morse oscillator which motivates 
our use of the McGehee transformation. 



For an excellent introduction to the BKL scenario we recommend a review paper 
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II. EQUATIONS OF MOTION 

A. Lagrangian 

An asymptotic form (near the cosmological singularity) of the dynamical equations of the 
non-diagonal Bianchi IX model reads (see, jl] and equations (2.25)-(2.28) in S*]): 

d'^lna b 2 d^lnb 2 ^ ^ 9^ In c 2 
dr'^ a ' dr'^ a b' dr'^ b 

The solutions to ([I]) must satisfy the condition: 

5 In a 9 In 6 9 In a 9 In c 9 In 6 9 In c 2 ^ ^ , , 

dr Ot dr dr Or dr a b 

It is easy to verify that ([1]) can be obtained from the Lagrangian equations of motion 

(ir((9x/) dxi' ' 

where Xi := In a, X2 := Info, X3 := Inc, and Xj := dxi/dr, and where the Lagrangian L has 
the form 

L := i;ii;2 + X1X3 + X2X3 + exp(2xi) + exp(x2 - Xi) + exp(x3 - X2). (4) 

B. Hamiltonian 

The momenta, pi := dL/dxj, are easily found to be 

p^=X2+X3, p2=Xi+Xs, p3 = X^+X2. (5) 

The Hamiltonian of the system has the form: 

1 . .1 

,Pi+P2+P3)~6xp(2xij-exp(x2-Xij-exp(a;3-; 

(6) 



H ■=piXi-L = -{piP2+PiP3+P2P3)-^{pl+pl+pl)-ew{'2xi)-exp{x2-Xi)-exp{x3-X2), 



which due to the relations 1^ and ([2]) leads to the dynamical constraint: 

H = 0. (7) 

The Hamilton equations 

Xj = {xj, H} = dH/dpi, pj = {pj, H} = -dH/dqj, (8) 

where 

/ d- d- d- d- 

y^\dxidpi dpidxj 



(9) 



3 



have the following explicit form: 



il = ^{-Pl+P2+P3), 

X2 = ^{Pl-P2+P3), 

is = ^(Pi +P2 -ps), 

Pi = 2 exp(2xi) — exp(x2 — xi), 

P2 = exp(x2 - xi) - exp(x3 - X2) 

P3 = exp(x3 - X2). 



(10) 

(11) 

(12) 

(13) 
(14) 
(15) 



Taking derivatives of ([T0l)- (fT2l) and making use of (fT3l) - (fT5l) leads directly to Eqs. ([I]), which 
proves that Lagrangian and Hamiltonian formulations are completely equivalent. 

The system (fT0l) - (fT5l) presents a set of nonlinear coupled differential equations. The 
space of solution of the above dynamical system is a 6D region of the phase space 
{xi,X2,X3,pi,p2,P3) G M^. This region is bounded by the constraint equation (JT]). Solv- 
ing ([7]) with respect to X3 gives 



X3 



X2 + log 

Substituting into (P|)-(ra we get: 



Pi P1P2 
4 2 



P2 



P1P3 ^ P2P3 



P3 



(16) 



Xi = 


-{-P1+P2+P3) 




X2 = 


^(Pl -P2 +P3), 




Pi = 


2g2xi _ ^-xi+X2 




P2 = 




P? 

4 


P3 = 


_g2xi _ ^-xi+X2 _ 


P? 
" 4 



P1P2 pl 

2 4 



P1P3 



P2P3 p| 
2^4' 



P1P2 _ P| , P1P3 P2P3 _ p| 
2 4 2 2 4' 



(17) 

(18) 
(19) 
(20) 

(21) 



Analytical solution to this system may not exist. For our purposes it is sufficient to 
understand the local properties of the phase space of the system by using the method of 
nonlinear dynamical systems (g-IO . 



III. DYNAMICAL SYSTEM ANALYSIS OF PHASE SPACE 

The local geometry of the phase space is characterized by the nature and position of 
its critical points^. These points are locations where the derivatives of all the dynamical 
variables vanish. These are the points where phase trajectories may start, end, intersect, 
etc. The trajectories can also begin or end at infinities. Then, after a suitable coordinate 
transformation projecting an unbounded phase space into a compact region (e.g. the so 



^ In what follows we use the terms critical or fixed points. 
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called Poincare projection), these originally infinite critical points become well defined for 
further analyses. The set of finite and infinite critical points and their characteristic, given 
by the properties of the Jacobian matrix of the linearized equations at those points, may 
provide a qualitative description of a given dynamical system. 

The above situation is specific to the case when a fixed point is of the hyperbolic type. 
In the case of the nonhyperbolic fixed point, linearized vector field at the fixed point cannot 
be used to characterize completely local properties of the phase space. 



A. Critical points of the unconstrained system 

Inserting Xi = = X2 = X3 = pi = p2 = ps into l.h.s. of (|TQ|) -f fT5|) leads to the following 
set of equations: 

= ^(-Pi+P2+P3), (22) 

= ^(Pi-P2+P3), (23) 

= l{pi+P2-P3), (24) 

= 2 exp(2xi) — exp(a;2 — a^i), (25) 

= exp(a;2 - a;i) - exp(a;3 - X2), (26) 

= exp(x3 - X2), (27) 

plus the Hamiltonian constraint equation: 

^{PiP2 + P1P3 + P2P3) - ^{pI+pI+pD -exp{2xi) - exp(x2 -Xi) -exp(x3 -X2) = 0. (28) 
Equations (I22])-(|24]): 

= +P2+P3, (29) 

= P1-P2 +P3, (30) 

= -P3, (31) 

(32) 

yield 

Pi = = p2= P3- (33) 
Equations fl25|) - fl27|) lead to the conditions: 

= exp(2xi), (34) 

= exp(x2 — xi), (35) 

= exp(x3 - X2), (36) 

Conditions dMD-dM]) are fulfilled for: 

xi — !■ —00 (37) 

X2 — Xi — i- —00 (38) 

X3 — X2 — )■ —00 (39) 
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thus for: 



Xi — i- — oo 

X2 — — OO, X2 < Xi < 0, 

X3 — )► —00, X'i < X2 < 0. 
Thus, the set of critical points fulfills the following conditions: 



Pi = = P2 =P3, 

xi —7- —00, X2 — 7- —00, X3 — 
X3 < X2 < xi < 0. 



-00, 



(40) 
(41) 
(42) 



(43) 
(44) 
(45) 



One may easily verify that this set satisfies the Hamiltonian constraint (!28|) . 
Thus the set of critical points Sb is given by 



{(a;i,a;2,a;3,pi,p2,P3) e M 
A(x3 <X2<Xi < 0); 



I (xi —J- —00, X2 — )■ —00, X3 — )■ 

=P2 =P3}, 



-00 



(46) 



where M := M U {-00, +00}. 

Stability properties are determined by the eigenvalues of the Jacobian of the system 
f|T0|) - f|T5|) . More precisely, one has to linearize equations f|T0|) - f|T5|) at each point. Inserting 
X = xq + 6x, where x = {xi,X2,X3,pi,p2,P3), and keeping terms up to 1st order in 6x leads 
to an evolution equation of the form 6x = J6x. Eigenvalues of J describe stability properties 
at the given point. 

The Jacobian J of the system f|T0|) - f|T5|) . evaluated at the set of critical points S, reads: 



J 



The characteristic polynomial associated with Jacobian J is equal to: 

p(A) = A^ 

and the eigenvalues are the following: 

(0,0,0,0,0,0). 



/O 


-1/2 


1/2 


1/2 \ 





1/2 


-1/2 


1/2 





1/2 


1/2 


-1/2 


























\0 








/ 



(47) 



(48) 



Since the real parts of all eigenvalues of the Jacobian are equal to zero, the set (H6|) 
consists of nonhyperbolic fixed points. 
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B. Critical points of the constrained system 

Inserting Xi = = X2 = pi = P2 = Ps into l.h.s. of f|T7|) -f l2T]) leads to the following set of 
equations: 

= ^(-Pi+P2+P3), (49) 

= ^(Pi-P2+P3), (50) 
= 2e2^'i -e-^■l+^^ (51) 
n _ ^2x1 , , pI P1P2 , pI PiPz P2P3 pI 

= -e^-i - e-^^+^^ _pi^pm_pi^pm^pm_pi_ ^^3^ 



Equations (pjl-ilSO]): 



= \{-pi+P2 + P3): (54) 
= \{pi-P2+p-i). (55) 



yield: 



Pi = P2, (56) 
P3 = 0. (57) 

Equations (!5T|) - (!53|) and the above solutions (!56l) - (l571) lead to the following set of equa- 
tions: 

= 2e2^i -e-^^+^^ (58) 
= e^^^ + 2e-^^+^^ (59) 
= -e^^'^ -e-^^+^^ (60) 

Solutions of the above set are the following: 

xi — )■ — oo (61) 
X2 — Xi — )■ — oo (62) 

(63) 

or equivalently: 

xi -> -oo (64) 

X2 — — OO, X2 < Xi. (65) 

(66) 

Thus, the set of critical points of the constrained system (fT7|) - (l2T|) reads: 

Pi =P2, P3 = 0, (67) 
xi — )■ — oo, X2 — — OO, (68) 
X2<xi< 0. (69) 
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Taking into account the Hamiltonian constraint ([6]) leads to additional condition 

X3- X2 ^ -00, (70) 

or, equivalently, to: 

X3 — )• —00, X3 < X2 < 0. (71) 

Requiring that the time derivative of 0:3, Eq f|T6|) . determined from the Hamiltonian 
constraint ([6]) satisfies f|T2l) , leads to Eq which together with fH9l) and fl50l) gives finally: 
Pi = = p2 = P3, instead of fl67j) . 

This way we have shown that in both cases of 'constrained' and 'unconstrained' dynamics 
the set of critical points S is defined by Eq. fHB]) . 

The Jacobian J of the system f[T7|) - fl^ evaluated at the set of critical points S reads: 



/o 


-1/2 


1/2 


1/2 \ 





1/2 


-1/2 


1/2 


























\0 








J 



The characteristic polynomial associated with Jacobian J is equal to: 

P(A) = A^ (72) 

so the eigenvalues are the following: 

(0,0,0,0,0) (73) 

which is consistent with the result obtained for the case of the unconstrained dynamics. 

In both cases of constraint and unconstraint dynamics we are dealing with the nonhy- 
perbolic type of critical points. Thus, getting insight into the structure of the space of 
orbits near such points requires an examination of the exact form of the vector field defining 
the phase space of our dynamical system. The information obtained from linearization is 
inconclusive. 

The results we have obtained so far have the following properties: 

1. The phase space is /izg'/ier dimensional. 

2. The critical points we are dealing with are of nonhyperbolic type. 

3. The set of critical points Sb, defined by fH6|) . is not a set of isolated points, but a 
3-dimensional continuous subspace of M^. 

4. The critical subspace 5*^ is situated in an asymptotic region of phase space with infinite 
values of its variables. 

Lower dimensional phase spaces can be easily analyzed by making mapping to known 
solved cases available in textbooks. Our six dimensional phase space needs sophisticated 
tools to be used. The nonhyperbolicity means that one cannot avoid direct examination of 
the original nonlinear set of equations defining the dynamics. The corresponding linearized 
set of equations is unable to reveal the nature of dynamics in the neighborhoods of fixed 
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points. The vector field may bifurcate there. However, suitable methods seem to be available 
(see, e.g. 0, IE!) to deal with such type of criticality. The third property is challenging. It 
needs discovering new methods within dynamical systems theory. Nonisolated critical points 
are atypical ones. The standard methods presented in textbooks [6-8] cannot cope with 
such a case. Making perturbations around critical points at infinities is not a well defined 
procedure. However, this problem may be avoided by mapping unbounded subspaces onto 
finite ones so the fourth property should not lead to serious problems. 



IV. THE MCGEHEE TYPE VARIABLES 

In this section, we apply the McGehee type transformations which has been success- 



fully applied to analyses of dynamical systems of general relativity (see, e.g. [12|, ll3|). In 
particular, we try to extend McGehee's transformations applied to the analysis of the Morse 
oscillator. The functional form of the vector field specific to Morse's oscillator looks similarly 
to our vector field f|T0|) - f|T5|) . This is why we hope to get profits using the McGehee mapping, 
which in the case of the Morse oscillator turns nonhyperbolic fixed points into hyperbolic 
ones. For completeness, we present in Appendix A the case of the Morse oscillator^. 
We define the transformation 

3 {Xi,X2,X3,Pi,P2,P3) > {Ui,U2,Us,Vi,V2,V3) G IR+ X (74) 

as follows: 

{ui,U2,U3,Vi,V2,V3) := (cxp Xi, exp X2, exp X3, pi, ps, Ps) • (75) 
We would like to interpret the dynamics in the new coordinates in terms of our original 



coordinates. This is possible if the mapping is canonical [141. Il5l|: 

{uk,ui}x,p = = {vk,vi}:^^p, {uk, vi}:^^p = 6ki, (76) 

where 

One may easily verify that 

{Uk, Ul}x,p = = {Vk, Vl}a:,p, {Ufc, Vl}x,p = Uk^kl- (78) 

Therefore, the transformation (ITlll is not canonical. Thus, the phase space in the new 
variables corresponds to a different dynamical system, i.e. the transformation maps the 
phase space of our original system to the phase space of a new system. As long as we stay 
at the classical level, i.e. we do not plan making quantization of our classical system, such 
procedure may be reasonable. It is so if the new system turns out to have simpler dynamics 
than the original one. In the rest of the present section we verify this expectation. 



For an introduction to this problem we recommend item 11 on page 722 of 
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A. Vector field 



Under this transformation, the vector field f|T0|) - f|T5|) turns into 

2ui = ui{-vi + V2 + vs), (79) 

2U2 = -V2+ V3), (80) 

2m3 = U3{vi +V2- fs), (81) 

Vl = 1u\ - U2/U1, (82) 

V2 = U2/U1 - U3/U2, (83) 

V3 = U3/U2. (84) 

To have well defined vector field, possibly of class , we rescale the time variable as follows: 

dr ds dr U1U2 ds U1U2 
Making use of ([85]) turns ([79])- ([84]) into: 

2u\ = ulu2{-vi + V2 + V3), (86) 

= uiul{vi - V2 + v^), (87) 

2^3 = UiU2U3{vi + V2 - Vs), (88) 

v[ = 2u\u2 - ul, (89) 

v'2 = u\ — U1U3, (90) 

= UlU-i. (91) 

We also redefine the constraint by multiplying it by the factor UiU2- The new form of 
the constraint reads: 

= U1U2H =: K, (92) 

where 

K = ^UiU2{ViV2 + V1V3 + V2V3) - ^MiU2(fi +vl+ fg) - ulu2 - u\ - U1U3. (93) 



B. Critical points 



All considerations carried out so far have been done with the assumption that Ui > 0,U2 > 
and U3 > 0. From now we lift this assumption down, i.e. we may have Mi > 0, ^2 > and 
M3 > 0. One can say that we extend the phase space of our system to a bigger one. 

Let us identify the critical points of the vector field (l86l)-(l9T]). We find that the set of the 
critical points is defined to be: 

Sm ■= {{Ul, U2, U3, Vl, V2, W3) G I Ml = = U2, M3 > 0, (f 1, V2, v^) E R^}. (94) 

It is clear that any element of Sm satisfies the constraint ( l93l) . One may verify that the 
transformation ( I7il) does not map the critical subspace Sb into Sm- 
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Let us examine the type of criticality specific to the set of critical points (IMj) . It is not 
difficult to find that the Jacobian of (l86l) - (!9T!) at any point of the critical set Sm reads: 

/ 0\ 





00000 
-Ms 
\ U3 00000/ 

The characteristic polynomial associated with J reads: -P(A) = A^. Thus, A = is an 
eigenvalue of J of multiplicity 6. Therefore, the set Sm consists of nonisolated nonhyperbolic 
fixed points. 



C. Hamiltonian structure 

Suppose that (I93p is a Hamiltonian of some new system. Let us derive the corresponding 
Hamilton's equations. One gets: 

(95) 
(96) 
(97) 
(98) 
(99) 
(100) 

which looks similarly to the system (15^ - fl9ip . but is still different. Thus, fl93l) cannot be the 
Hamiltonian of the original vector field f l86p - fl91l) . 



2uiu[ 


= Ulu2{-Vi +V2+ V3), 


2U2U'2 


= UiuKvi -V2 + V3), 


2M3% 


= UiU2U3{Vi +V2 -fs) 


uiv[ 


= 2ulu2 — 


U2V2 


= ul- U1U3, 


M3M3 


= U1U3, 



V. THE POINCARE TYPE VARIABLES 

Since examination of phase space at 'infinite region', f H6|) . is difficult mathematically, 
we change coordinates of the phase space to map the set of critical points ( H6l) onto a 
finite region. We map the infinite space into a finite Poincare sphere, parameterized by 
Cartesian coordinates {Xi, X2, X3, Pi, P2, P3), as follows: 



Xi =: 


Xi 

1 — r ' 


(101) 


X2 =: 


X2 

1 — r ' 


(102) 


X3 =: 


X3 

1 — r ' 


(103) 


Pi =■■ 


Pi 
1 — r ' 


(104) 


P2 =■■ 


P2 
1 — r ' 


(105) 


P3 =■■ 


P3 
1 — r ' 


(106) 
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where = Xf + X| + X| + + P2 + ^3^. We also rescale the time parameter r by defining 
the new time parameter T such that dT := (ir/(l — r). In these coordinates our phase space 
is contained within a sphere of radius one - 'infinities' correspond to r = 1. 
If the mapping is canonical, we should have: 



(107) 



{Xl, Xk}x,p — — {Pi, Pk}x,p, {Xl, Pk}x,p 

The map (I10ip - fll06p is not canonical, because we have: 

{Xk,Xi}x,p = (1 - r)g{a){xkPi - Xipk), (108) 

{Pk, Pi}x,p = (1 - r)f{a){xkPi - XiPk), (109) 

{Xk, Pi}x,p = (1 - rYhi - (1 - r) {f{a)xkXi + g{a)pkPi) , (110) 

where a := r^/(l — r)^, /(a) 7^ 0, g{a) ^ 0. It is clear that there is no chance to satisfy 
f ll07p for any r including the limit r — )■ 1. 

The insertion of f lTIIT]) - ffTUB]) into ([10])- ([15]) gives: 



1 — r 

X2 
1 — r 

^3 

1 — r 

Pi 
1 — r 

1 — r 

1 — r 



-Pi + P2 + P3 



1 



P2 
P2 



P3), 



1^1 — r) (2exp Y 
1^1 — t) ( exp 
1^1 — r) (exp 



2X1 



exp 



X2 — Xi 



X2 — Xi 



exp 



1 — r 
X3 — X2 - 
1 — r ' 



X3 — X2 
1 — r 



(111) 
(112) 
(113) 
(114) 
(115) 
(116) 



where prime denotes derivative with respect to the new time parameter T. 

To find the fixed points we insert X[ = = X^ = X^ = P^ = P^ = P^ into ffTTT|) - ffTTB|) 
by using the elementary formulas: 



-^r = {X,X[ + X^X'^ + X3X^ + P,P[ + P2P2 + P^Pz)!^ 



and, e.g. 



d ^ X, \ _ X[{1 
dT\l-r) ( 
After rearrangement of terms we finally get: 



xy 



(117) 
(118) 



PI + P2 + P3 
P1-P2 + P3 
P1 + P2- P3 



2 exp 



2X1 



exp 



1 — r 
X2 — Xi 

1 — r 



exp 
exp 
exp 



X2 




1 


— r 


^3 


-X2 


1 


— r 


X3 


-X2 



0, 


(119) 


0, 


(120) 


0, 


(121) 


0, 


(122) 


0, 


(123) 


0. 


(124) 
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The solution to f lTT9|) - fn:2T|l reads: Pi = = Pa = Ps- The equations ffT22D- ffT2il) can be 
satisfied in the hmit r i— )■ 1 if 

'2Xi X2 — Xi ^3 — X2 
nm exp = = hm exp = hm exp , (125) 

r^-l- 1 — r r-s-1- 1 — r r->l- 1 — T 

which leads to the condition: X3 < X2 < Xi < 0. Therefore, the critical subspace is defined 
to be: 

Sp := {(Xi,X2,X3,Pi,P2,P3) I (X3 < X2 < Xi < 0) A (Pi = = P2 = P3)}. (126) 

It is not difficult to verify that the transformation f ll0ip -( |T06l) does not map Sb into Sp. 

It is clear that any point of Sp, in the limit r — )■ 1~, satisfies the constraint ([7]) which in 
the variables fll0ip - fll06p has the form: 

^^^(PiP^ + P1P3 + - (A^ + Pi + Pi) (127) 

2X1 X2 — Xi X3 — X2 

— exp exp exp = 0. (128) 

1 — r 1 — r 1 — r 

One can resolve (either manually or by symbolic computations) the nonlinear vector 
field flllip - flll6p with respect to the derivatives X(, Xg, . . . , P3, and find the corresponding 
Jacobian. Its value at any point of the subspace Sp (in the limit r 1— )■ 1) turns out to be 
a six dimensional zero matrix. It means that linearization of the exact vector field, at the 
set of critical points Sp, cannot help in the understanding of the mathematical structure of 
the space of orbits of considered vector field. An examination of the nonlinearity cannot be 
avoided. One may say, formally, that the set Sp consists of the nonhyperbolic type of fixed 
points. 



VI. THE DIRAC TYPE VARIABLES 



Solutions to Eqs. (niI]l - (fT^ define the kinematical phase space J-^. The physical phase 
space J-'p is defined by those solutions to (fTUi) - (IT5ll . which additionally satisfy Eq. ([7]). 

The physical Dirac observable O is defined to be a function on the phase space that weakly 
Poisson commutes with first-class constraints of the dynamical system [isl ^ 20| . In our case 



there is only one constraint ([7]) so the equation defining the physical Dirac observables reads 
{0,H)^,, where {,.}:. X: (|-|- " |-|-). im 



k=l 



and where the sign f» means that (11290 has to be treated as a weak equation in Dirac's 
sense Functions which satisfy f ll29p strongly define the kinematical observables. So 

they are defined by the solutions to the following equation 

[-pi +P2 +P3)^ h (pi -P2+Pd,)-^ y- {P1+P2 -p?,)^— 

0x1 0x2 ox^ 

8(0 do do 

+2(2e^^i - e^^-^^) — + 2(e^^-^i - e"^)— + 26^^-^^ — 
opi dp2 ops 

= 0. (130) 
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The physical observable can be obtained from the solutions to Eqs. (11301) satisfying the 
constraint ( 128|) . 

For the purpose of quantization, which we plan to carry out in the future, we rewrite the 
classical dynamics in terms of the Dirac observables. More precisely, we give the physical 



meaning to the following factorization, consistent with the Darboux theorem [15|, |16 



(^dxk A dpk^ = (^dqa A diTa^ — dt A dH, (131) 

k=l a=l 

where {pm, x„) G J-'k, which can be made in the neighborhood of the constraint surface H = 0. 
It is clear that the transformation {xi,X2,X3,pi,p2,P3) — > {qi,q2,t,ni,7i2, H) underlying 
(I13ip is canonical, but not unique. 

Now, suppose that [qa^i^p) are the Dirac observables. On the constraint surface, due to 
(I13ip and considerations on time- dependent vector fields and differential forms in ISHlTJ, 
we propose the following: 

d d dH 

= =: {g„, i^Tjg.TT - -j^qa ^ -i?^qa = {qa, HT}g,n = -J. — (132) 



and 



d d OHt 

= {tTo, -ffjl^^o =: {tTc, HT}g,n - -J^T^a =^ -7:7; = {iXa, Hrjq^T, = j: , (133) 



where 



^ / d- d- d- d- \ 



a=l 



- dqa dTTa Oil a dqa 



Therefore, the dynamics is generated by the Hamiltonian Hrp (to be found) and is pa- 
rameterized by an evolution parameter T. In this new setting the system has no dynamical 
constraints. One may examine the above dynamics using the dynamical systems methods 
0,9. The phase space is now only four dimensional, which simplifies analysis of the original 
dynamics, f|T0l) - f[T5l) . defined in six dimensional phase space. 

The factorization procedure described above make sense, but not in the neighborhood 
of the subspace of the critical points Sb defined by fH6|) . and it is not defined uniquely. 
However, we expect that taking into account the possible structure of the set of critical 
points of the vector field (I132p - (ll33p may emerge a natural choice for Ht and T. 



VII. CONCLUSIONS 

The novelty of our approach is atypical formulation of the Hamiltonian description of the 
non-diagonal Bianchi IX dynamics. Our point of departure is taking the Lagrange equations 
in the asymptotic region near the cosmo logical singularity. Next, we make an educated guess 
what the Lagrangian and Hamiltonian could be which would lead to these equations by using 
the principle of least action. Usually one starts with an exact form of the Lagrangian (or 
Hamiltonian) and analyzing dynamics one gets an asymptotic form of the dynamics. Our 
approach guarantees that we are really dealing with the prototype of the BKL scenario as 
the BKL theory was derived by generalization of the Bianchi IX model dynamics considered 
in this paper [1, 0]- 
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Taking into account the Hamiltonian constraint in two ways, before or after identifying 
the critical points, leads to the same results concerning the set of critical points of the vector 
field of the considered dynamical system. 

The circumstance that critical points occur in asymptotic regions seems to be not an 
obstacle because one may project an unbounded subspace onto a finite region. 

The original phase space presented in Sec. Ill is /iz^f/ier dimensional. Since a// eigenvalues 
of the Jacobian (corresponding to the nonlinear vector field) are purely imaginary, no reduc- 
tion to lower dimensional phase space is possible, because the center manifold theory 0, Isj 
cannot be applied. Since all the critical points are nonhyperbolic, the information obtained 
from linearization is inconclusive. One has to apply, e.g. the theory of normal forms 0, [1] 
to get insight into the structure of the space of orbits in the neighborhoods of the critical 
points. 

Additional difficulty results from the fact that the set of critical points is not a set of 
isolated points, but a 3-dimensional continuous subspace of R^. Perturbation of the vector 
field in the neighborhood of such a space, seems to require making analyses in transverse 
space to the space of critical points. 

The mapping of the original system by making use of sophisticated transformations into 
other systems (helpful in analyses of other gravitational systems) does not simplify the 
problem. The critical points are still nonhyperbolic. This is the case of both the McGehee 
and the Poincare transformations. In the former case, the 'miracle' of getting away from the 
nonhyperbolicity observed in the case of Morse's oscillator does not occur. Apart from this, 
these mappings are obtained by non-canonical transformations. Such type of transformations 
of phase space lead to different dynamical systems so quantization of the dynamics before 
and after such transformation would give different results. 

We believe that parameterizing the phase space by the Dirac observables, proposed in Sec. 
VI, is the best choice. Since the physical phase space defined there is only four dimensional, 
we expect that the subspace of critical points has even lower dimensionality so it can be 
examined by the phase portrait methods 0, TO]- Finding the Hamiltonian that generates 
the dynamics which is no longer the dynamical constraint, outlined in Sec. VI, may enable 
the construction of the quantum evolution of the Bianchi IX universe. Preliminary results 
obtained in the papers 17, 21-2^ for the FRW and the Bianchi I models which make use, 
to some extent, the above ideas are encouraging. 

The realization of the factorization described in Sec. VI and quantization of the resulting 
theory may be done, but it is beyond the scope of the present paper. This intriguing issue 
will be the subject of our next paper. 
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Appendix A: Morse oscillator 

In this section we consider the Morse oscillator^. It is defined by the Hamiltonian: 

H{x,y) = y^-P{e--le-'n, (Al) 
where {x,y) G M? and < /3 G M is a constant. The Hamilton equations read: 

X = dH/dy = y, (A2) 
y = -dH/dx = -(3 (e"" - e'^-), (A3) 

where dot denotes derivative with respect to the time parameter t. The critical points of 
the vector field (lA2|) - (lA3p . obtained by solving x = = y, are easily found to be {xi,yi) = 
(0,0) =: Ci and {x2,y2) = (oo,0) =: C2. 

The characteristic equation of the Jacobian at ci, corresponding to (lA2|) - (]A3p . reads: 



= det{j{ci) - XI) = det ^ }^ ^ 



(A^-/3), (A4) 



and has the solution 

(-V^,V^), (A5) 

which means that the critical point at ci is of a hyperbolic type since the real parts of both 
eigenvalues are different from zero. 

The characteristic polynomial corresponding to the critical point at C2 has the form: 



= det{j{c2) - AI) = det ^ }^ ^ 



A^ (A6) 



with the solution 

(0,0). (A7) 

Since both eigenvalues are purely imaginary, the critical point at C2 is of a nonhyperbolic 
type. 

Now, let us apply the McGehee type transformation of the variables |8|: 

u := e-"/^ (A8) 
V := y. (A9) 

and rescale the time as follows: 

dt/ds = -2/u. (AlO) 
The equations (]A2|) -( lA3|) are mapped into the system: 

u = V, (All) 
V = 2l3{u-u^), (A12) 

where prime denotes derivative with respect to the new time parameter s. 



It is a special case of the so called driven Morse oscillator 
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It is clear that the vector field (lAlip - (lA12l has three fixed points: {ui,Vi) = (0,0) =: 
fi, {u2,V2) = (0,1) =: /2 and (M3,f3) = (0,-1) =: /s. The Jacobian corresponding to 
(lATT]) -(lAT2l reads: 

J{u,v) = I ^ OA i ) • (A13) 



2/3(1-^2) 

It is easy to verify that J{fi) = J{f2) = Jifs)- The characteristic polynomial equation, e.g. 
for J{fi), reads: — 2(3 = 0, and has the solution 

^,VW). (A14) 



Thus, all three critical points are of hyperbolic type. 

It is interesting to examine the role of the rescaling of the time flAlOl) . Without the 
rescaling, the mapping flA8l) - flA9l) turns the system flA2p - flA3P into the vector field: 

ii = -uv/2, (A15) 
V = (3(u^-u^)^ (A16) 

which has three critical points at: (0,f) =: ei, (—1,0) =: 62, (1,0) =: 63. The first one is 
in fact a manifold as G M. The Jacobian corresponding to the vector field f lA15|) - (lA16P 
reads: 

The characteristic polynomial det(^J{ei) — Al) = has two solutions (0, — f/2), so ei is a 
nonhyperbolic fixed point. Since J{e2) = —J{es), the characteristic polynomials in both 
cases have the same roots, which read: {—i^/]3,i^/]3). Thus, both 62 and 63 are centers. We 
conclude that the time rescaling is crucial. 
One may easily verify that we have: 

{u, u}^^y = = {v, v}^^y, {u, v}^^y = -u/2, (AlS) 

which means that the transformation (1A8I) - (1A9I) is not canonical. 

Let us verify if the new vector field (IAlll) -( lA12l) can be given the structure of the Hamil- 
tonian vector field. It is not difficult to find that (lAlip - (]A12p can by obtained via Hamilton's 
principle from the Hamiltonian K{u, v) defined to be: 

K{u, v) := v'12 - P{u^ - u^/2) = /2 - /3(e" - 6-^72) = H{x, y). (A19) 

Thus, the transformation (IA8l) - flA10P leads from one Hamiltonian system to another one, 
and the relation between both systems reads: 

yx — H{x, y) = vu — K{u, v). (A20) 

The relation f lA20p looks like a standard relation between two Hamiltonian systems related 
by canonical transformations with the generating function equal to zero p^], but now the 
time derivatives on both sides of (1A20P are determined with respect to different variables 
related by (lATOj) . 
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